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ABSTRACT 



Diagnostic model output parameters, provided by the Fleet 
Numerical Weather Central, Monterey, California (FNWC) , and 
the marine fog frequency climatology developed at the Naval 
Postgraduate School, Monterey, California, are statistically 
processed in context with marine surface synoptic ship re- 
ports for the purpose of developing a linear regression 
scheme for modeling the distribution of marine fog. The study 
area encompasses a large section of the North Pacific Ocean 
(from 30-60N) . The time period involves 0000 GMT, 1-30 July 
1976. The predictand is a fog parameter developed from a 
quantitative categorization of each of the 4481 synoptic ob- 
servations according to reported present and past weather, 
low-cloud type, and visibility codes. The 38 model output 
and climatological predictor parameters are interpolated to 
each of the ship observation positions and the resultant data 
file is used for the derivation of the regression equations. 
The diagnostic capabilities of the regression equations, 
along with other existent approaches, are analyzed through the 
use of three verification scoring systems--Heidke Skill, 
Threat, and the Panof sky-Brier Probability scores. Improve- 
ment over climatology and FNWC's operational fog probability 
program (FTER, applied at analysis time) is demonstrated. 
Selective mappings of the regression equation outputs and 
categorized observations are intercompared with the sea-level 
pressure analysis; FTER; and the most significant predictor 
parameter in the regression equations, evaporative heat flux. 
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I. INTRODUCTION AND BACKGROUND 



Marine fog continues to be a threat to safe nautical and 
low-level aviation activities, civilian and military. Wheeler 
and Leipper (1974) documented various monetary and human 
losses associated with United States Navy operations in poor 
visibilities due solely to fog. An accurate analysis of 
marine fog, and its associated increases in forecasting 
accuracy, should help decrease these losses as well as en- 
hance the specifications of low-level cloud inputs to the 
various planetary boundary layer models in the environmental 
sciences community. 

Research on marine fog by the Departments of Meteorology 
and Oceanography at the Naval Postgraduate School (NPS) has 
concentrated on the climatology and the diagnosis of this 
phenomenon from observational data (Renard, Englebretson , and 
Daughenbaugh , 1975; Renard, 1976). A recent effort (Ihli and 
Renard, 1977) attempted to utilize satellite data in the 
diagnosis of marine fog regimes. The results indicate that 
the approach has potential, but more data are needed for fur- 
ther testing. In the realm of forecasting marine fog, an 
earlier study (Nelson, 1972) attempted an approach similar to 
the Model Output Statistics (MOS) method developed by the 
National Weather Service (NWS) (Glahn and Lowry, 1972) . 

Nelson correlated visibility at sea over the North Atlantic 
Ocean to locally observed parameter values and developed 
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regression equations using the "perfect prognostic" approach. 
Out of the 20 or so parameters tested, wind speed, relative 
moisture content, and parameters associated with evapora- 
tion showed the highest correlations with visibility. He then 
tried to forecast visibility using, as predictors, grid point 
output values of these parameters from the operational models 
at the United States Navy Fleet Numerical Weather Central, 
Monterey, California (FNWC) . The results showed little skill 
in the forecast mode. 

Only FNWC is presently producing fog forecasts on a hemi- 
spheric basis through their statistical probabilistic product 
called FTER (see Appendix A) . This product is the result of 
calculations based on a statistical processing of certain 
model output parameters related to fog occurrence (U. S. Naval 
Weather Service Command, 1975) . FTER does not use any type 
of climatological parameter, partly due to the unavailability 
of accurate climatologies at the time of the development of 
the product. Only a limited evaluation of the FTER product 
has been accomplished to date (Renard, 1975b) . 
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II. OBJECTIVES AND APPROACH 



The primary objective of this study is to utilize the MOS 
concept in the development of a statistical equation (or set 
of equations) that would accurately specify the distribution 
of marine fog over an extensive oceanic area, such as the 
North Pacific Ocean. The secondary objective is to evaluate 
the skill of an existing operational scheme, FNWC ' s FTER 
product, and the NPS fog frequency climatology for the same 
region (Englebretson , 1974; Daughenbaugh, 1975; Willms, 1975). 
In essence, this is an "imperfect diagnostic" approach de- 
veloped in part from numerically analyzed/predicted atmos- 
pheric and oceanographic parameters, similar to the MOS 
concept. 

The approach used in pursuit of the primary objective 
involved three steps: a) the classification of each synoptic 

observation in the data base (described in Section IV) by 
cloud, weather, and visibility reported (predictand) ; b) the 
interpolation of model output grid-point parameters (and 
other parameters derived therefrom) and a climatological 
parameter (predictors) to the observation points; and c) the 
search for statistical relations between these predictors 
and the predictand. The immediate goal was to find one or 
more multiple, linear regression equations that would be 
significantly more skillful than either the NPS climatology 
or the FNWC FTER product in describing the actual distribu- 
tion of fog at analysis time. 
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III. DATA 



A. AREA OF STUDY 

The NPS climatological fog frequency study (Willms, 1975) 
was confined to the North Pacific Ocean. Since this was 
judged to be the most accurate climatology of fog occurrence 
to date, the same general area was chosen for this study. 
Further restrictions were dictated by the desire to avoid 
the influences on fog patterns by large land masses and the 
large-scale upwelling in the waters off of the western coast- 
line of the United States. The area of study is displayed in 
Fig. 1 on a Mercator projection. The area is outlined using 
the boundary I,J points of the standard 63-by-63 grid used 
by FNWC. Note the resultant curvature of the grid lines when 
the grid is placed on a Mercator projection. Fig. 2 has the 
same area and outline displayed on a Polar Stereographic pro- 
jection. The entire FNWC I,J grid, on a Polar Stereographic 
projection, is shown in Fig. 3 for reference. 

B. TIME PERIOD FOR STUDY 

The month of July, 1976 was chosen as the basic time 
period. The main reason for this selection is the high 
relative frequency of marine fog for this month, as displayed 
in the NPS climatology (Willms, 1975). 

The 0000 GMT synoptic time was chosen as best since this 
time is local noon over the international dateline (180 de- 
grees) with daylight hours occurring over all of the area of 
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study. This factor should give maximum relative accuracy of 
observations by the transient ship observers (Nelson, 1972). 

Other aspects involved in using only 0000 GMT data are 
persistence and diurnal variations. It is widely held that 
diurnal variations of meteorological phenomena in oceanic 
regions are minimal. In any case, the use of one synoptic 
time will eliminate any diurnal effect. 

Panofsky and Brier (1968) suggest that, in order to 
eliminate persistence at stationary observation points, data 
from every third or fourth day be used in the dependent 
sample. Since the vast majority of the observations used for 
this study x^ere from transient (non-stationary) ships, such 
persistence effects are considered negligible. 

C. SYNOPTIC WEATHER REPORTS 

FNWC provided all of the primary- time synoptic weather 
observations from their data archives. Duplicate reports and 
reports obviously in error were eliminated from the observa- 
tions provided. In all, 4481 0000 GMT observations were ob- 
tained from the period 1-30 July 1976. Observations for 
31 July 1976 were not available from FNWC. 

While the majority of the observations came from tran- 
sient ships, Ocean-Weather-Station P and 11 land (island or 
coastal) stations (see Fig. 1) were also included in the 4481 
observations. The 11 land stations are listed in Appendix B, 
each with its station number, elevation, and location (lati- 
tude, longitude). The location of land station number 32618 
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was adjusted from its actual location (see Fig. 1) to just 
inside the study area. 

D. MODEL OUTPUT PREDICTOR PARAMETERS 

FNUC also provided all of the selected model output para- 
meters (MOP's) for the time periods involved. These MOP's, 
their description and/or applicable remarks, and range of 
values (from histograms not shown) are listed in Appendix C. 
Some of the MOP's were chosen after consultation with Mr. Leo 
Clarke, Director of Research, FNWC, and others were chosen 
using the experience and advice of earlier studies (Schramm, 
1966; Nelson, 1972; Grisham, 1973, Renard 1975a). The only 
desired MOP not available from the FNWC archives is a dew- 
point variable, but the EAIR (vapor pressure) parameter used 
in the study is derived from the dew point. The acronyms 
or abbreviations listed in Appendix C will be used as the 
variable names or symbols throughout the remainder of the 
report . 

E. DERIVED PREDICTOR PARAMETERS 

Other parameters used in the study were derived from the 
basic set of 21 MOP's provided by FNWC.''' These are listed 
in Appendix D in the same manner as the basic MOP set is 
listed in Appendix C. As with the basic set, the acronyms 
and abbreviations listed are used in the remainder of the 
report . 

"'"Except Sensible Heat Flux (SHF) which was recovered 
from SHF = SEHF - EHF ; it is listed in Appendix C. 
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F. CLIMATOLOGICAL FOG FREQUENCY PREDICTOR 

The values of the NPS July climatological fog frequencies 
(Willms, 1975) used in the study were determined from inter- 
polations to a 249-by-249 grid, which is about 1/4 of the 
FNWC I,J grid (Fig. 4). The resolution afforded by this 
fine spacing enabled the researcher to define all of the major 
features on the NPS July climatology charts. Such spacing is 
compatible with the original development of the fog fre- 
quencies, which were organized by one-degree-by-one-degree 
grid squares. Some subjective extensions to Willms' (1975) 
final isopleths were necessary in order to accommodate to 
the study area. These extensions are shown as the dot-dash 
isopleths in Fig. 4. 

The overlay used for determining the climatological values 
at the grid points of the 249-by-249 grid was developed by 
computer plotting latitude, longitude values from a modified 
FNWC I, J-to-latitude/longitude program (See Appendix E) . 

The other predictor variables were available for the standard 
63-by-63 FNWC I,J grid-point locations. The abbreviation, 
CLIMO, used in the following portions of the text, stands for 
the NPS fog frequency climatology. 
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IV. FOGCAT PROGRAM 



A. BACKGROUND 

Until a methodology for using satellite data to identify 
marine fog regimes becomes operational, the main information 
sources for this purpose are synoptic ship reports and some 
island station reports. Previous researchers at the NPS used 
three elements of a ship's synoptic weather observations to 
categorize or differentiate between reports indicating fog 
("foggers") and those indicating weather other than fog ("non- 
foggers") . The three elements used were present weather (ww) , 
past weather (W) , and visibility (VIS), See Tables I and II. 
While these elements certainly can be used to ascertain the 
occurrence/non-occurrence of fog, at least one additional 
element appears useful in determining the existence of fog 
regimes, namely, low cloud type (CL). See Table I. 

Some of the combinations of ww, W, and CL used here as 
indicators of the occurrence of marine fog generally would 
not be considered valid in most verification schemes. How- 
ever, this non-verification is primarily due to the existing 
rules in weather observing (U. S. Department of Commerce, 

1969) . These rules require the reporting of precipitation 
(ww more than 49, such as light drizzle or rain), even though 
a co-existing non-precipitation phenomenon (such as one of 
the modes of fog in the ww 40s group) may have more operational 
significance . 
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B. CATEGORIZATION PROGRAM 



A scheme (FOGCAT; Fog Cat egorization) was devised to place 
each of the 4481 observations into one of the five major cate- 
gories: Strong foggers (S) , Foggers (F) , Past/weak foggers 

(P) , Maybe foggers (M) , or Non-foggers (N) . The scheme 
evaluates only the present weather (ww) , past weather (W) , 
and low cloud type (CL) codes for category placement. The 
synoptic visibility (VIS) code is not considered until the 
score assignment phase. A description of the FOGCAT scheme 
follows . 

The 72 combinations of ww, W, and CL that are used in 
categorizing the observations are listed in Table IV, using 
the legend listed in Table III. The partitioning of the 
major categories into subcategories was done for three pur- 
poses, the least important of which is the collection of 
occurrence statistics for each combination. The listing order 
of the subcategories in Table IV is a subjective estimate of 
the degree of certainty of a given combination indicating the 
presence of operationally important fog at an observation 
point. The listing order was arranged with a view to both 
the skill of the observer and the precedence rules for report- 
ing. In effect, it is assumed that the more and/or the 
denser the fog indicators reported by an observer (trained or 
untrained) in an observation, the more certain it is that the 
observation was taken in or near an existent fog regime. The 
second aspect of the listing order is that it provides a 
means of assigning a probabilistic score to each such 
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observation of ww, W, and CL; that is, a numerical value of 
the degree of certainty that fog was occurring. 

The main problem encountered in devising FOGCAT involved 
the grouping of the various values of ww, W, and CL. The 
groupings used are listed in Table III, which is also the 
legend for Table IV. There are several major differences 
between the present (ww) and past (W) weather groupings 
therein and those used in the development of the NPS fog fre- 
quency climatology (Willms, 1975). The 50's (drizzle) and 
some of the 60 's (rain) code values (see Table I) for ww are 
separated from the "any others" group. The ww code values 
of 10 (deep light fog) , 28 (heavy fog in past hour) , and 40 
(light fog) are placed in the same group (10G) . The code 
values of 20 (drizzle in past hour) and 24 (freezing drizzle 
in past hour) for ww are moved from the "any others" group to 
the one (11G) containing the values of 11 and 12 (shallow 
heavy fog) . The elements in each of the latter two groupings 
do not necessarily have similarities in time or nature of fog 
but are related herein primarily due to a subjective consid- 
eration of their respective indicative qualities for the 
diagnosing of fog regimes. 

The W code value of 5 (drizzle) is paired with 4 (fog) 
and this group (4,5 in Table III) is treated as a reliable 
fog indicator only when in combination with the ww groups of 
41G (deep heavy fog) or 10G (the latter only for major cate- 
gory S) . Otherwise, the W value of 5 is grouped (5G) with W 
values of 0-2 (sky coverage changes) . 
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The low cloud (CL) type groupings are more straightfor- 
ward. Values 5 (s tr a to cumulus) and 7 (fractos tratus of bad 
weather) are used to counter the possibility of an untrained 
observer substituting either of these for the more properly 
reported 6 (stratus or fractostratus) . CL values reported 
as (/) indicate an obscuration. Of the 4481 observations 
used in the study, 919 have a CL of (/), and 466 of these are 
in the a2 subcategory. 

1 . Major Category S 

The ww, W, and CL combinations comprising the strong 
fogger category generally have at least two strong fog indi- 
cators in combination. The strongest combination indicating 
fog is given by subcategory al. However, over 50 per cent of 
the total S cases, 793, are categorized in the a2 subcate- 
gory. A relatively high number of all foggers were expected 
to fall into this subcategory because of the obscuring nature 
of fog and the results of Willms ' (1975) study. 

2 . Major Category F 

In general, there are at most two fog indicators in 
the combinations comprising this category. Placing subcate- 
gory dl into the F, rather than S, category indicates the 
author's opinion that the CL value of B is not as certain a 
fog indicator as the CL value of 6 in subcategory cl. Sub- 
categories d3 and d4, which are the first occurrences in the 
listing order of a 50G ww value, can be considered equivalent 
to al and a2 in a fog situation except that the precedence 
reporting rule prohibits reporting one of the 41G values. 
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But the total number of such cases is small and the effect 
on the final regression equation should also be minor. 

3 . Major Category P 

The past/weak fogger category consists mainly of com- 
binations which have a W value of 4 and no strong indication 
of fog at or within the past hour of observation time. 

4. Major Category M 

The "maybe" category represents the gray area of the 
scheme. Some of the reports that are placed in this category 
are likely to be in or near a fog regime, but are not veri- 
fiable as fog. 

5 . Major Category N 

The non-fogger category was expected to include the 
majority of all the cases, and it did, as evidenced by 2604 
observations or 56 per cent of the total. 

C. SCORE ASSIGNMENT PHASE 

In order to develop a multiple linear regression equation 
for specifying the probability of fog occurrence, it was 
necessary to assign a quasi-continuous numerical value to 
each unique combination of ww, W, and CL in the synoptic re- 
port. The nature of the FOGCAT table facilitates this value 
assignment, as well as the assignment of discrete, numerical 
or yes/no designations. 

The continuous value assignments to the FOGCAT elements 
are illustrated in Table V. The letter portion of each sub- 
category symbol in Table IV was used to position the ww, W, 
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and CL combinations into the listing order of Table V. For 
example, the four subcategories, bl thru b4, are combined into 
"b". This grouping of the four combinations that comprise 
subcategories bl thru b4 effectively gives each of these 
four combinations the same degree of certainty. 

The subcategory letter is paired with the visibility (VIS) 
value of each observation and a corresponding Base Score is 
determined. For example, an observation that has a VIS value 
of 92, and which receives a subcategory placement of kl or k2, 
would be assigned a Base Score of 360. As such, Table V was 
used as a decision logic table. 

It can be seen that a particular Base Score is assigned 
to many unique combinations of observed values of ww, W, CL, 
and VIS. For example, the Base Score of 260 can be assigned 
to any of 15 such combinations (five subcategories times the 
three visibility groups) . The non-usage of the Base Scores 
of 160 and 180 makes a distinct gap or buffer zone in the 
continuous value range between those scores assigned to non- 
foggers and the ones assigned to observations with any 
possibility of being a fogger. 

The assigned Base Score was then scaled by dividing it 
by 5.8. This yielded the Fog Score, which was used as the 
predictand (dependent variable) during the regression phase. 
The Fog Score values (rounded) are listed next to the Base 
Score values in Table V. 
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D. FURTHER REMARKS ON FOGCAT 

Other researchers would undoubtedly use different place- 
ments of the combinations of ww, W, and CL in the listing 
order of Table IV. Still others might suggest that the CL 
values could be used only if the sky cover code (N) value is 
co-evaluated. The Base Score assignment procedure could be 
kept "pure" by having a Base Score for each of the 72 combi- 
nations listed in Table IV and restructured in Table V. The 
Base Score range could then be 0-770 (intervals of 10 rather 
than 20), and/or the intervals could be variable so that Base 
Scores for foggers are close together but those for non- 
foggers would be farther apart (given consecutively listed 
subcategories). But, as a means to accomplish what is 
basically a feasibility study, FOGCAT is felt to be a "reason- 
able" scheme. 

Support for Nelson's (1972) general comments on the 
accuracy of the transient ship synoptic observations is shown 
in the counts listed in Table IV. According to the observing 
rules (U. S. Department of Commerce, 1969), those observations 
having present weather values in the 40 ' s should also have VIS 
values in the 90-93 range. Subcategory a2 ended up with a 
total case count of 466. By the rules, all 466 of these ob- 
servations should have ended up being listed in the "Poor VIS" 
count for the a2 combination. Only 360 were so listed. It 
would have been possible to "bogus" in the proper visibility 
for those 106 faulty observations, but this was not attempted. 



The number of observations in Table IV for the a2 sub- 
category and the N major category have already been discussed. 
At least one other such count deserves mention, namely, the 
count of 313 observations for subcategory nl. The meaning 
of this unexpected high count is unclear, but the CL value 
of 6 suggests these observations may have been taken at or 
near fog regimes. A shift in the position in the listing 
order of Table IV for this subcategory should definitely have 
a noticeable effect on the resultant regression equations. 
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V. PREDICTOR INTERPOLATION 



The interpolation scheme chosen to provide values of the 
38 predictors (21 MOP's, 16 derived MOP's, and CLIMO) to the 
observation points is a natural bicubic spline curvilinear 
interpolation scheme, locally called SPLIN. The locally de- 
veloped computer program is available at the NPS W. R. Church 
Computer Center. The scheme can handle large arrays (up to 
300-by-300) . But for time considerations, it was decided to 
limit the array to a 4-by-4 matrix for each interpolation 
performed. 

In order to use any interpolation scheme, it is necessary 
to be able to find the observation points in the same frame 
of reference as the predictors. This was accomplished through 
the latitude/longi tude- to-I , J- coordinates conversion formulae 
described in Appendix F. The resultant I,J coordinates 
(63-by-63 grid) for a given observation point were used to 
generate the coordinates required with respect to the 249 -by - 
249 grid used to derive the CLIMO values . 

The nature of the interpolation scheme, which is essen- 
tially a cubic curve fitting process, is such that errors can 
be introduced when the majority of the 4-by-4 array values are 
very similar to one another in sign and magnitude. These 
errors, and their possible effects, will be discussed later. 
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VI. REGRESSION SCHEME 



A. BACKGROUND 

A stepwise multiple linear regression program, called 
BMDP2R (University of California, 1973), is available at the 
NPS W. R. Church Computer Center, and was chosen as the means 
to derive a marine fog diagnostic equation. BMDP2R computes 
a sequence of multiple linear equations in a stepwise manner. 
At each step one variable is added to or deleted from the 
previous step's equation as dependent on the F-to-enter and 
F-to-remove criteria. Generally, the variable entered into 
the equation, in any given step, is the one which makes the 

greatest addition to the amount of predictand variance ex- 

2 

plained (AR ) . 

The BMDP2R program permits the specification as to wheth- 
er or not the y-intercept (YI) is to be treated as a variable 
or predictor. It is also possible to specify the tolerance 
level, which is a factor controlling the possible linear in- 
terdependencies or correlations between the various predictor 
variables. This tolerance specification is used by the pro- 
gram to inhibit the entering of a variable which has a high 
correlation to another previously entered variable. A case 
in point may be seen in Table VI, which lists correlations of 
all predictors to Fog Score along with the ranges of the 
variables. It can be seen that the correlations between the 
Fog Score predictand and both the EHF and SEHF predictors are 



similar in magnitude and sign. The first regression equation 
obtained (see Table VII) used EHF as a predictor, but not 
SEHF. The correlation coefficient between EHF and SEHF is 
0.94 for the data used, which shows the high, and expected 
interdependence between these two predictor variables, and 
thence, why the tolerance control will tend to keep SEHF out 
of the equation. The other correlations between the predic- 
tor variables are not shown. 

The approach taken in this study was to allow the treat- 
ment of the y-intercept (YI) as a variable and to use the 
BMDP2R program's default values of 4.0, 3.9, and 0.01 as the 
values for the minimum acceptable F-to-enter, maximum accept- 
able F-to-remove, and the minimum acceptable tolerance 
criteria, respectively. 

B. OUTPUT RESULTS 
1 . Correlations 

The simple correlation coefficients (cc) , between Fog 
Score and the 38 predictor variables, for the entire set of 
4481 observations are listed in Table VI. The groupings of 
the variables are for easier comparison between variables 
with similar characteristics, such as temperature, and/or 
model source. See Appendices C and D for model sources. 

The largest absolute cc value is the one listed for 
the EHF predictor. All of the other heat-flux type variables, 
except SOLARAD, also have relatively high cc values. 

Also worth noting are the ccs for the FTR and CLIMO 
predictor variables. This leads to the expectation 
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that the results in Section VII (Verification Scoring) would 
show that FTER possesses little, if any, skill over CLIMO in 
the diagnosis of the July 1976 marine fog regimes. 

2 . Regression Equations 

The FTER and DDWW variables were removed from the 
BMDP2R program's use list before the generation of the first 
set of regression equations was developed. FTER was removed 
because it is the existent predictor over which improvement 
is desired. DDWW was removed because it is a cyclic (0 = 360 
degrees) variable. Cycle variables are ill-conditioned for 
use in regression programs. Other cyclic variables, which 
were not recognized as such until late in the study, are PDW 
and SDW. Their inclusion as variables in some of the resul- 
tant equations may be undesirable, but time restrictions pro- 
hibited regeneration of the equations without them. 

The set of equations generated by BMDP2R for the 

whole study area (4481 cases) is listed in Table VII. They 

are listed in their respective stepwise order along with the 

variable regression coefficients, YI values, explanations of 
2 2 

variance (R ), and the increase in R at the end of each step. 
A discussion of Eq. (9.9) will occur later in the verification 
section. 

It is readily seen that the EHF variable, in conjunc- 
tion with YI, is the most significant predictor of the Fog 
Score. Since SEHF and the rest of the heat flux parameters 
were excluded by the tolerance specification, a separate run 
was made with EHF removed from the program’s use list, with 
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the purpose of determining whether SEHF would be entered 
first in place of EHF. This was expected, because SEHF has 
the second largest cc magnitude. The results (not shown) 
upheld the expectations. SEHF did replace EHF as the first 
variable to enter. 

Along with the EHF variable, only CLIMO, EAIR, and 

TSEA are significant variables in the explanation of the Fog 

2 

Score variance. This significance is based on an R increase 
of 0.005 rather than the widely held standard of 0.01. If 
0.01 were used as the significance limiter, EHF would be the 
only such variable, along with the constant YI. 

None of the variables listed for Eq. (9) were unex- 
pected. All of them can be logically and physically associa- 
ted with the occurrence of marine fog. The association of 
T925 and fog is rather difficult to understand. This corre- 
lation is probably primarily due to the meridional distribu- 
tions of fog regimes and air temperatures, as is known to be 
the case with TSEA. 

Two other equations generated for the whole area are 

listed in Table VIII. The final step equation and its respec 

2 

tive coefficients are listed, as well as the AR due to each 

2 

variable entered, the final R , variables deleted (if any) 

2 

and their associated AR at entry and removal. Variables not 
given consideration (not in use list) during the regression 
are also listed. 

The development of Eq . (9T) stems from a transforma- 
tion of the Fog Score variable (predictand) . This 
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transformation (T) takes the form of Fog Score = Fog Score - 

CLIMO. Effectively, this forces the regression coefficient 

of CLIMO to be one (1,0), as is reflected in Table VIII. 

Thus , CLIMO was not allowed entry into the equation as a pre- 

o 

dictor during the stepwise procedure. The final R" for Eq. 

(9T) indicates that this approach holds little promise for 
useful application. 

The second equation listed in Table VIII is Eq. (9NC) . 
In this experiment CLIMO was removed as a predictor variable. 
This regression result was run at the request of FNWC due to 
unavailability of acceptable fog frequency climatology for 
all oceans. The results, except for entry order, are similar 
to Eq. (9). The PDW variable is higher on the selection list 
than its counterpart in Eq. (9) and its significance is 
higher. In any case, the influence of this cyclic variable 
is minimal. 

Two sub-area equations were also obtained. The whole 

area file was sorted and copied into two smaller files, one 

[for Eq . (39)] for those cases (1504 observations) south of 

40N, and the other [for Eq. (40)] for all the other cases 

(2977 observations), namely, at or poleward of 40N. The final- 

step forms of these two equations are listed in Table IX. The 
o 

final R for Eq. (40) was the highest obtained, while that for 
Eq. (39) was the lowest, with the exception of Eq . (9T) . This 
"poor” showing by Eq. (39) might be attributed to the low 
occurrence of fog in the southern zonal belt below 40N. There 
were only 107 foggers (major category S and F) in the 1504 
observation reports in that area. 
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Two other runs made with BMDP2R are also worth men- 
tioning. One run had FTER as the only predictor variable in 

the use list, the other run had CLIMO as such. The runs were 

2 

made against the whole area data file. The final R for the 
FTER run is 0.5092, and for the CLIMO run, 0.5391. This 
strengthened the expectation that FTER would possess little, 
if any, skill over CLIMO as a fog predictor. 
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VII. VERIFICATION SCORING 



A. BACKGROUND 

Two types of verification scores are used in the study to 
test the "skill" of different predictors in describing the 
distribution of the FOGCAT major category designators (S, F, 

P, M, N) . The first type of scoring is used in judging the 
effectiveness in forecasting ("nowcasting" in this study) 
discrete (event, non-event) occurrences of a given phenomenon. 
The second type of scoring is used to judge the effectiveness 
of probabilities in estimating the occurrence of the event. 

In all, three scores are used in this study. The first two, 
Heidke Skill Score (HSS) and Threat Score (TS) , are of the 
first type of scoring approach. The third score, Probability 
Score (PS) , evaluates the probability of fog as assigned by 
Fog Score, FTER, or CLIMO, and is of the second type described 
above . 

The exact formulae used for the three scores (HSS, TS , and 
PS) are given in Appendix G. The Heidke Skill Score has been 
used for quite some time in the meteorological community but 
may be falling into disfavor currently, while the Threat Score 
is currently in vogue with the National Weather Service. The 
Probability Score is appropriate to a predictor given in terms 
of a probability with a range of 0-100 per cent. 
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B. RESULTS 

1 . Heidke Skill and Threat Scores 

Tables X through XIII show the various results for the 
HSS and TS calculations. The FOGCAT major categories of S 
and F were used as the criteria for classifying an observa- 
tion as a "fogger". The cutoff or cut value for a predictor 
is considered to be that value at which the estimate of the 
event is changed from a "yes event" estimate to a "no event" 
estimate (or vice versa) . Only those cut values listed for 
FTER and CLIMO were used, because of the results of an 
earlier study (Renard, 1975b) . The cut values for the equa- 
tions are those giving the best verification results. 

Table X reflects the results hopefully expected at 
the beginning of the study, namely, Equation (9)'s relatively 
significant improvement over either FTER or CLIMO as a marine 
fog diagnostic tool. CLIMO and FTER are quite comparable 

in skill, which upholds the earlier expectation. Eq . (9T) , 

2 

which shows "poor" R results, is comparable to Eq. (9) in 
HSS and TS results, as are Eqs. (9NC) and (9.9). The listing 
of the scores for Eqs. (1) through (8) shows the generally 
step-by-step increase in the scores. This listing complies 
with a suggestion by Panofsky and Brier (1968) to test each 
of the equations in the stepwise regression scheme. 

Listed in Table XI are the number of observations per 
day as well as the associated skill scores for the three pre- 
dictors FTER, CLIMO, and Eq. (9), using the best cutoff 
values from Table X for each of them. The day-by-day 
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comparison shows a good improvement of Eq.(9) over the capa- 
bilities of FTER and CLIMO in the diagnosing of the event. 

The day-by-day comparisons in Table XII for Threat Scores 
further highlights the merits of Eq. (9) in comparison to 
FTER and CLIMO. 

The results listed in Table XIII for the two zonal 
belts again show relatively significant improvement by the 
regression equations over FTER and CLIMO. The equations 
are comparable within the respective areas. That is, Eq. (9) 
(whole area equation) exhibits about as much skill as the 
equations generated specifically for the sub-areas. 

2 . Probability Scores 

These scores for the three predictor types are listed 
in Tables XIV and XV. CLIMO shows a slight improvement over 
FTER (the smaller the value, the better the skill) while 
Eqs . (9), (9T) , and (9NC) show very little improvement over 
CLIMO. Eq. (9.9) shows the best improvement. This equation 
was obtained by changing the YI value of Eq. (9) with the 
express goal of optimizing the resultant PS. However, the 
optimization effect did not carry over to the HSS and TS 
results listed in Table X. 

The zonal-belt PS results listed in Table XIV show 
that the equations have a slight improvement in the northern 
belt over FTER and CLIMO. But in the southern belt CLIMO has 
the best score. This is probably due to the preponderance of 
"good" weather (non-fogger) observations south of 40N. 
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The daily PS values listed in Table XV support the 
findings in Table XIV, namely, that Eq. (9) does not show 
much improvement over CLIMO and/or FTER. The daily perform- 
ance of Eq. (9.9) (optimized) shows a far better than average 
improvement over Eq. (9) , FTER, and CLIMO on a temporal fre- 
quency basis. 

The relatively small improvement of Eq. (9) over CLIMO 
and FTER in the PS results, as compared to the TS and HSS 
results, can be examined more closely by using Table XVI. 

Eq. (9) does a better job of specifying the foggers in the 
higher probability classes than CLIMO or FTER while CLIMO 
and FTER perform better in indicating non-foggers in the 
lower probability classes. The results of the PS calculations 
indicate these two forecasting aspects yield comparable PS 
values for CLIMO, FTER, and the regression equations. The 
figures listed in Table XVI compare favorably to those of 
Renard (1975b) for July 1974 data. 

Whenever a predictor probability value was negative, 
zero was substituted so as not to cause complications with 
the PS formula. (See Table XIV.) This adjustment was not 
necessary for the HSS and TS calculations. 
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VIII. COMMENTS 



In section V it was indicated that the interpolation 
scheme might introduce errors into the data base. This error 
production can be readily seen by comparing the value ranges 
of the predictor parameters before and after the interpola- 
tion. The "before" values are listed in Appendices C and D, 
while the "after" values are listed in Table VI. The effect 
of this error introduction is clearly shown by the cases of 
below-zero value for the FTER variable listed in Table XIV. 
The FTER range has been expanded from 0.0 to -0.09 on the low 
end and from 1.00 to 1.02 on the high end. The values may 
not seem significant, but 1183 observations were assigned 
negative values of FTER by the interpolation process rather 
than the more desirable zero or small positive values. To 
check the extent of this effect, two additional Probability 
Score runs were made in which each of the 4481 values of 
FTER was modified. The results of these runs are the FTER 
+0.02 and FTER +0.05 entries in Table XIV. It can be seen 
that 891 of the 1183 observations received interpolated FTER 
values of 0.0 to -0.02, which is judged to be an acceptable 
deviation from zero. 

Mappings of sea-level pressure, EHF, FTER, probabilities 
from Eq. (9), and verifying observations by FOGCAT are dis- 
played for days 12 and 24 in Figures 5 through 14. The EHF 
charts are included since this model output parameter is 
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predominant in explaining a large amount of the variance of 
fog. Day 12 was chosen because the scores produced for FTER 
and Eq. (9) are about "average", considering all days 
studied, while the performances of these two predictors for 
day 24 are "well above average". 
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IX. CONCLUSIONS AND RECOMMENDATIONS 



The primary objective of this study was to utilize the 
Model Output Statistics approach in the development of a sta- 
tistical regression equation that can accurately specify the 
distribution of marine fog. Through the interpretation, 
categorization, and statistical processing of various data 
elements, an equation was obtained whose variables are 
logically and/or physically related to fog (Table VII). The 
secondary objective, that of evaluating the skill of the 
different predictors (FTER, CLIMO, and the equations obtained) 
against a given data base, was fully achieved. 

The methodology of using diagnostic model output para- 
meters, rather than observational values, for describing the 
initial distribution of marine fog has tremendous potential. 
The equations obtained show significant increases in skill in 
performing this task in comparison to the use of climatology 
or the current FNWC fog forecasting product, FTER. FTER was 
shown to have no appreciable skill over the NPS fog frequency 
climatology (CLIMO) for the criteria used. 

Regression Eq. (9), developed by using the whole area 
data file, proved to be as skillful in zonal band application 
as were the equations generated for the zonal bands . This 
suggests that the zonal band partitioning does not segregate 
dissimilar areas well enough with respect to fog occurrence. 
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The following recommendations are offered for both short 
range and long range continuations of this work: 

(1) Process the July 1976 parameters and observations 
again, using one or more different interpolation 
schemes . 

(2) Break up the area of study into a multiple number of 
sub-areas on both geographical/meteorological bases 
and derive regression equations for each area. 

(3) Recover the long wave (LW) portion of the THF para- 
meter and redo the July 1976 investigation including 
this parameter. 

(4) Test the regression scheme in a forecast mode, using 
primitive equation prognostic model output parameters. 

(5) Perform the same procedures, and modifications listed 
in these recommendations, on data from different 
months and/or synoptic times. 

(6) Process the data bases using a discriminate analysis 
procedure . 

(7) Analyze scatter diagrams from the data base using a 
non-linear approach. 

(8) Plot the FOGCAT designations for each synoptic 
report each day and research the areas of fog 
generation, maintenance, and advection. Relate 
these fields to changes in the patterns of signifi- 
cant variables such as EHF. 

(9) Experiment with changes to the FOGCAT structure and 
Fog Score assignment method or intervals. 
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( 10 ) 



Use Eq. (9) or (9NC) on a test operational basis as 
soon as practicable, and test the skill against FTER. 

(11) Introduce a fog persistence parameter into any re- 
gression scheme for prognostic probabilities. This 
parameter also may be useful for the diagnostic 
probability equations. 

(12) Eliminate the advection terms from consideration 
during the regression runs and redo the July 1976 
investigation. This would decrease the complexity 
of the parameters required for the "final" product. 
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APPENDIX A - Description of FNWC 1 s Fog Probability 
Forecast Program (U. S. Naval Weather 
Service Command, 1975) . 

1 . Program Development 

The forecast for Advection Fog probability is based 
on a system of multi-parameter tests, where each parameter 
is weighted on the basis of test results. The weights of 
all parameters are summed, divided by 5, and normalized to 
give a value of zero to one, representing probability. 0.1 
would be 10% probability; 0.9 would be 907<>. 

Advection Fog has five tests: 

(1) When surface air temperatures are near 0°C, 
there is a high probability of fog. (Based on S. Petterssen) 
Weights favoring fog are assigned in the + 5°C range from 
0°C. 

(2) When condensation is taking place, there is 
usually an effective fog condition. This test is done by 
using evaporation calculations from moisture parameters. If 
evaporation is negative or zero, a certain weight is assigned 
to allow for the condensation factor. 

(3) If the near surface relative humidity exceeds a 
certain critical value (> 95%) , a high weight is assigned on 
fog probability. 

(4) An air mass discriminant is applied. If the 
surface geostrophic wind is being advected toward warmer 
(or colder) surface conditions, the air mass would be 
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classified accordingly- -as in the old k and w air-mass system. 
A warm air mass would call for a high weight. The weight 
assigned is based upon calculations of the air mass wind com- 
ponent normal to the isotherms. 

(5) A weight is derived from application of a Fog 
Criterion. The criterion is based upon the change in dew 
point that occurs during the mixing of two air masses of dis- 
similar character. The resulting dew point depression in the 
mixed air is related to fog probability, which establishes a 
certain criterion value, a strong probability weight is 
assigned. 

User activities receive the fog probability forecasts 
in the form of a field, where each grid point has a fog 
probability value. 

2 . Program Input 

(1) Sea-surface temperature analysis. 

(2) Surface air temperature analysis. 

(3) Surface vapor pressure analysis. 

(4) Sea-level pressure 6-hour ly analysis. 

(5) P.E. prognoses of the above parameters from 0-48 

hours . 

3 . Program Computation 

The various tests are calculated as described above 
for parameters at Tau's equal to 0, 6, 12 to 72 hours. The 
fog probability FTER is calculated for all land and sea 
points . 
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4. Program Limitations 

The five-parameter test method is admittedly some- 
what gross in character, but has enough built-in controls 
to insure a correct answer--given correct humidity values. 

The program gives excellent results in cases of broad-scale 
fog. The weakness of the program stems almost entirely 
from use of the large 63x63 grid. Local effects just cannot 
be accounted for, and fog is particularly susceptible to 
small-scale effects such as cold air drainage, small bodies 
of water, and man-made air pollution. Forecasters, in using 
the Advection Fog products, should always make allowances for 
known local effects. 
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APPENDIX B: Land Stations Used in Study (U. S. Air Force, 

1972) . 



ion Number 


Latitude 


Longitude 


Surface Ele' 
(ft) 


32174 


45.2 


N 


147.9 


E 


38 


32186 


46.2 


N 


150.5 


E 


73 


32195 


46.9 


N 


151.9 


E 


26 


32207 


48.3 


N 


153.3 


E 


55 


32213 


50.9 


N 


156.7 


E 


42 


32215* 


50.7 


N 


156.2 


E 


unknown 


32217 


50.0 


N 


155.4 


E 


11 


32559 


53.1 


N 


160.0 


E 


88 


32594 


51.5 


N 

■k-k 


156.5 


E 

'k’k 


6 


32618 


54.9 


N 


166.2 


E 


19 


70454 


51.9 


N 


176.7 


W 


14 



Station not listed in reference. Lat/long values listed 
are as given with each synoptic report furnished by FNWC. 



"'Location given slightly removed from actual location of 
55.2 N latitude, 166.0 E longitude to accommodate study 
area. 
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APPENDIX C. Output Parameter Descriptions. 



Source Model 



Symbol Name 
Description 



: k k k 

Mean Minimum Maximum Units 



A. Analysis Parameters (FNWC's Mass Structure Model) 



B. 



PS 


Sea Level Pressure: 


1017.2 981.2 


1034.0 


(mb) 




Analysis of observed 


sea level parameter. 




TAIR 


Surface Air 
Temperature 


17.9 3.5 


30.3 


(°C) 




Analysis of observation-level air temperature . 




EAIR 


Surface Vapor 
Pressure : 


18.1 3.0 


38.4 


(mb) 




Analysis of observation- level vapor pressure derived 




from the dew point. 








T925 


925 mb Air 
Temperature : 


16.2 -4.4 


26.4 


(°c) 




Analysis of 925 mb air temperature . 






TSEA 


Sea Surface 
Temperature: 


16.6 4.1 


29.0 


(°C) 




Once- daily analysis of observed sea- 


■surface temperature. 


P.E. 


Parameters (FNWC's Primitive Equation Model) (Kaitala, 


1974) 


TX 


Surface Air 
Temperature : 


18.5 5.3 


28.2 


(°C) 



EX 



Derived from surface air and potential temperatures, 
boundary layer depth, upper- level winds extrapolated to 
surface, air density, drag coefficient, gustiness factor, 
and empirical constants. 

Surface Vapor 



Pressure : 


18.4 


6.6 


33.7 


(mb) 


Derived from model ' s 


mixing 


ratio . 






SOLARAD Solar Radiation: 


46.5 


25.9 


88.4 


(g*°c / 










cm 2 *H) 



Calculated absorption of incoming short-wave (solar) 
radiation (positive downward) (H = depth unit) . 



EHF 



Evaporative Heat 
Flux: 



5.5 



■14.9 



26.3 



(g;°c/ 

cm 2 *H) 

Derived using air density, drag coefficient, extrapolated 
winds, and mixing ratios. 
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APPENDIX C (continued) 



SEHF Sensible/Evaporative 

Heat Flux: 1.9 -39.7 27.9 (g-°C/ 

cm 2 *H) 

SEHF = Sensible Heat Flux (SHF) + EHF 

SHF Sensible Heat Flux: -5.6 -28.4 3.6 (g*°C / 

cm 2 *H) 

Recovered from SHF = SEHF-EHF. Originally derived by FNWC 
using drag coefficient, extrapolated winds, surface air 
temperature, TX, density, and constants. 

THF Total Heat Flux: -37.6 -80.9 -9.0 (g*°C/ 

cm 2 *H) 

THF = SEHF - SOLARAD + LW, where LW is the heating due 
to long-wave (terrestrial) radiation. 

C. Marine Wind Model (FNWC) 



WWW 


Marine Wind Speed: 


13.3 


0.0 


74.8 


(knots) 


DEWW 


Marine Wind 
Direction : 


18.7 


0.0 


36.0 


(degrees/ 



10 ) 

Both variables derived from a dynamic balancing of surface 
wind and sea- level pressure. 

D. Spectral Ocean Wave t-fodel (S.O.W.M.) (FNWC) 



HW 


Significant Wave 
Height: 


4.4 


0.0 


22.1 


(feet) 


PPW 


Primary Wave Period: 


8.2 


0.0 


20.0 


(sec) 


PEW 


Primary Wave 
Direction : 


18.2 


0.0 


33.7 


(degrees/ 

10) 


SPW 


Secondary Wave 
Period: 


5.6 


0.0 


25.0 


(sec) 


SEW 


Secondary Wave 
Direction 


11.5 


0.0 


33.7 


(degrees/ 

10) 


WCP 


Probability of 
White Caps: 


0.3 


0.0 


7.0 


(per 

cent) 



APPENDIX C (continued) 

E. Other Model Output Parameters (FNlvC) 

SSTA Sea Surface Tempera- 
ture Anomaly: -0.6 -4.5 3.0 (°C) 

Calculated anomaly of sea- surface temperature from the 
mean of the day as interpolated from the monthly mean 
values . 

FIER Probability of Fog: 0.1 0.0 1.0 (per cent) 

Forecast of advecticn fog probability. 

* 

These values are for the grid points in the rectangular area (Fig. 3) 
encompassing the area of study for 1-30 July 1976, including 0.0 for 
wave parameters over land points. 
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APPENDIX D: Derived Parameter Descriptions 



Symbol 


Name 

Description 


k 

Mean 


k 

Minimum 


k 

Maximum 


Units 


u 


Zonal Wind 
Component 

u = -WWW sin (DDWW • 


1.7 

10). 


-20.9 


28.9 


(m/sec) 


V 


Meridional Wind 
Component 

v = -WWW cos (DEWW • 


1.3 

10). 


-25.2 


25.4 


(m/ sec) 


CAPU 


I Directional 
Wind Component: 

CAPU = -u • sin(LNGA) 


1.2 

- v • 


-24.8 

cos (LNGA), 


20.5 

(Haltiner, 


(m/ sec) 
1971) . 


CAPV 


J Directional 
Wind Component: 

CAPV = u • cos (LNGA) 
where LNGA = -10 -(I, 


-2.0 -35.1 

- v • sin (LNGA) , 
J point longitude) 


18.7 

(Haltiner , 


(m/sec) 
1971) . 


THETAX 


Potential 
Temperature X: 
Derived using PS, TX. 


290.1 


277.5 


300.7 


(°K) 


THETAR 


Potential 
Temperature R: 


289.5 


276.5 


301.7 


(°K) 



STABX 



STABR 



Derived using PS, TAIR. 



Stability X: -0.06 -0.19 0.07 

Derived using [THETAX -(THETA of T925) ]/ (PS-925) . 

Value greater than zero indicates absolute instability. 



(°K/mb) 



Stability R: -0.07 -0.25 0.12 

Derived using [THETAR -(THETA of T925) ]/ (PS-925) . 
Same value effect as STABX. 



(°K/mb) 



ASTDX 


Air- Sea Temperature 








ASTDR 


Difference X: 

ASTDX = TX - TSEA. 

i 

Air- Sea Temperature 


1.9 -3.4 


11.8 


(°C) 




Difference R: 

ASTDR = TAIR - TSEA. 


1.3 -5.6 


9.7 


(°C) 


ADTSEA 


Advecticn of TSEA: 

Formulae 


)_t , 

/W\A 

0.04 -0.79 

and notes below. 


0.95 


(°c/ 

Hour) 


ADTX 


Advection of TX: 


0.35 “''-0.73 


1.02 


(°c/ 

Hour) 



Formulae and notes below. 
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APPENDIX D (continued) 



ADTAIR 


Advection of TAIR: 0.04 

Formulae and notes below. 

JLJLJU 

a/CX 


-0.85 


0.86 


(°c/ 

Hour) 


AASTDX 


Advection of ASTDX: -0.00 

Formulae and notes below. 


-0.69 


0.78 


(°c/ 

Hour) 


AASTDR 


_ i _ t — * - 

/Owf 

Advection of ASTDR: 0.00 

Formulae and notes below. 


-1.02 


0.64 


(°c/ 

Hour) 



"These values are for the grid points in the rectangular area (see Fig. 3) 
encompassing the area of study. 

I,J point longitude calculated using I, J-to-latitude/longitude 
conversion formulae described in Appendix E. 



Mean may be biased towards zero. Only values used in study were 
calculated. Grid points outside study area received zero. 

Advection Formulae and Conditions: 

For the advection of a quantity (R) the formula, ADQ = -V*V(Q), was 
used in the following finite difference form: 

AK, = - ™»[CAFU- <Q I+1 - Qj.p j + CAFV- (Q J+1 - Qj.^] , 

where HMAP = (1 + sin(60))/(l 4- sin (latitude)) 

and EM = [ (2) • (6.37-10 6 ) • (1 + sin (60))]/31.205 

(31.205 = grid mesh lengths, pole to equator, on FNWC's I,J grid). 

In the temperature advection calculation for point C, using the five 
grid points illustrated below, one or two of the points, namely, A, B, 

D, or E, may be outside the study area. In the bogusing method 
suggested by Mr. Leo Clarke, FNWC, when a non-center point (e.g. , point A) 
was judged to probably produce a land and/or dissimilar sea area 
influence on the resulting advection, the center point value (point C) 
was substituted for it. This "bogusing" is necessary to maintain a 
"purely" marine characteristic in the resultant parameter value. 

The set of study area boundary grid points used for bogusing (and 
sane double bogusing) are depicted in Fig. 2. The upper case letters 
denote the points whose values were used for the bogusing. The lower 
case letters mark the positions of the adjacent points being bogused. 

The only boundary point close enough to land to give concern is to 
the one near station number 32594 (south tip , Kamchatka Peninsula) . 
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APPENDIX D (continued) 

However, it was judged to be over the ocean after exaimrung TX, TAIR, 
and TSEA values at this and adjacent land/ocean I,J-grid points. 

A 

• • O 

I+1,J 

BCD 
I,j+1 I*J I,J-1 

I 

• • • 

I-1,J 
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APPENDIX E: Conversion of I,J grid coordinates to latitude/ 

longitude positional values (FNWC Subroutine 
Library) . 

The following formulae were programmed in Fortran as a 
subroutine used by several different large programs developed 
for this study: 

LONGITUDE = K - arctan [ (J-J ) / (I-I ) ] 

(RED) 2 - (I-I D ) 2 - (J-J„) 2 

LATITUDE = arcs in 7 E-y 

(RED) Z + (I-I p ) Z + (J+J p ) 

where: K = -10 or 170 depending on quadrant used in 

study, (-longitude means West longitude) 

*.u 

J = Pole position = 31 
P 

it 

Ip = Pole position = 31 

RED = Distance from pole to equator in grid mesh 
lengths = 31.205 

*32 used on IBM 360/70 at NPS W. R. Church 
Computer Center. 
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APPENDIX F: Conversion of latitude/ longitude position to I,J 

grid coordinates (FNWC Subroutine Library) „ 



The following expressions were programmed for use during 
the interpolation phase of the study: 



I - I p + [RED • ( T p • cos (-10+X) ] 

J = J p + [RED • • sin(-10+X)] 

where = Latitude 

A = Longitude (negative = West) 

I = Pole position = 31 

P 

Vc 

Jp = Pole position = 31 

RED = Distance from pole to equator in grid mesh 
lengths = 31.205 



32 used on IBM 360/70 at NPS W. R. Church 
Computer Center. 
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APPENDIX G: Verification Score Formulae 



From a contingency table and associated information the 
Heidke Skill Score (HSS) and Threat Score (TS) can be 
calculated: 



EVENT ESTIMATED 



EVENT 

OBSERVED 



where 



YES NO 



YES A 

NO B 



C Total (T) =A+B+C+D 

No. of Correct 
D Forecast (FC) = A + D 



HSS 



FC - EX 
T - EX 



Range : 



-2BC 

2 

B z + C 



j £ HSS 



< 



1 



EX = 



(A+B) (A+C) + (D+B) (D+C) 

T 

correct forecasts due to 



No. of expected 
chance . 



TS 



A _ A 
T-D A+B+C 



Range : 0 £ TS £ 1 



Both scores indicate more skill with larger 
positive values. 



The Probability Score (PS) is from that given by 
Panofsky and Brier (1958) and may be written as (R.enard, 1975) : 

PS - £.£ Pi 2 ♦ E 1 a - Pj) 2 ) 

i-i j=i 

Range: 0 £ PS £ 2 

where N = Total number of cases 

n = Total number of non-events 
o 

= Associated probability value for the non-event 
n^ = Total number of events 

Pj = Associated probability value for the event 
The closer to zero, the greater the skill. 
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TABLE I. Abridged version of internationally used weather-code figures 
and definitions for reporting present and past weather and 
low clouds in the surface synoptic report (U. S. Departments 
of Comnerce, Defense, and Transportation, 1969). 



Present Weather 
Code 



Value 


Definition 


00-03 


Characteristic change 
of the state of the sky 
(cloud) during the past 
hour. 


04-09 


Haze, dust, sand, or 
smoke. 


10 


Deep light fog. 


11-12 


Shallow heavy fog. 


13-17 


lightning, thunder, or 
precipitation within 
sight, not reaching the 
ground. 


18-19 


Squall ( s ) , funnel 
cloud(s) during the 
past hour. 


20 


Drizzle during the 
past hour. 


21-23 


Rain, snow, or rain and 
snow during the past 
hour. 


24 


Freezing drizzle during 
the past hour. 


25-27 


Shower (s) during the 
preceding hour. 


28 


Fog during the past 
hour. 


29 


Thunderstorm during the 
past hour. 



Present Weather 
Code 



Value 


Definition 


30-39 


Dusts torn, sandstorm, 
drifting or blowing snow. 


40 


Fog at distance, but not 
at station during the 
past hour (visibility 
less than 1 km) . 


41-49 


Deep heavy fog at the 
time of observation 
(visibility less than 
1 kn) . 


50-59 


Drizzle, or drizzle and 
rain. 


60-63 


Slight to moderate rain. 


64-65 


Heavy rain. 


66 


Slight freezing rain. 


67 


Moderate or heavy 
freezing rain. 


68 


Slight rain or drizzle 
and snow. 


69 


Moderate or heavy rain 
or drizzle and snow. 


70-79 


Solid precipitation not 
in showers. 


80-99 


Showery precipitation 
or precipitation with 
current or recent 
thunderstorms . 



55 



TABLE I. (continued) 



Past Weather 


Code 


Value 


Definition 


0 


Cloud covering % or less 
of sky throughout the 
period. 


1 


Cloud covering more than 
% of sky during part of 
period. 


2 


Cloud covering more than 
% of sky throughout the 
period. 


3 


Sandstorm, or duststorm, 
or blowing snow. 


4 


Heavy fog, thick haze, 
or smoke. 


5 


Drizzle 


6 


Rain 


7 


Snow, rain and snow 
mixed, or ice pellets. 


8 


Shower (s) 


9 


Thunderstorm, with or 
without precipitation. 


Visibility 


90 


Less than 50 m 


91 


0-199 m 


92 


200-499 m 


93 


500 m - 0.99 km 


94 


1 - 1.99 km 


95 


2 - 3.99 km 


96 


4 - 9.99 km 


97-99 


equal to or greater than 
10 km. 



Low Cloud Type 

Code 

Value Definition 

0 No low clouds 

1 Ragged cumulus of fair 
weather. 

2 Generally towering 
cumulus. 

3 Cumulonimbus without 
cirriform or anvil tops. 

4 Stratocumulus formed by 
cumulus spreading out. 

5 Stratocumulus not formed 
by cumulus spreading. 

6 Stratus or fractostratus. 

7 Fractostratus of bad 
weather. 

8 Cumulus and stratocumulus, 
with bases at different 
levels . 

9 Cumulonimbus with 
cirriform top. 

/ Low cloud obscured. 
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TABLE II. Listing and description of the elements of the visibility-weather group in the 

synoptic reports used to identify the existence of marine fog in earlier studies. 
(After Wil lms, 1975). 
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q 
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d 
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Q 
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Q O 




o 


o 
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o o 
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•iM 


•H 


Mti *H 




•iM 


•H 


O 


0) 


•H *iM 
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4-1 


4-J 


q 


4-J 


4-1 


•H 


5m q 
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a> 


f 
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Cti 


CL) Cti 
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Horizontal W 90-94** < 2 

visibility 95 , 96** 2-10 

* At time of visibility restriction indicated by present or past weather code figure. 
'^"Visibility codes used to identify fog when ww(W) is not 10, 11, 12, 28, or 40-49(4). 



TABLE III. Groupings and symbols used in FOGCAT 
categorization scheme. See Table IV. 



Present 


Weather (ww) 


Symbol 


Associated 
ww Codes 


41G 


41-49 


10G 


10,28,40 


11G 


11,12,20,24 


50G 


50-59 


60G 


60-63,66,68 


* 


any ww not 
listed above 



Past 


Weather (W) 


Symbol 


Associated 
W Codes 


4 


4 


4,5 


4,5 


5G 


0,1, 2, 5 


2G 


0,1,2 


JL> 

A 


any W not 
listed above 



Low Cloud (CL) 

Symbol Associated 
CL Codes 

6 6 

5,7 5,7 

B / 

* any CL not 

listed above 



G = Group. B = Low clouds obscurred. 
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TABLE IV. 



Scheme for categorizing observations according to 
likelihood of fog (FOGCAT) . See Table III for 
symbols. Number of July 1976 North Pacific Ocean 
observations in each subcategory stratified by 
visibility code. 



Major 


Sub- 


Present 


Past 


Low 

Cloud 


Number of Visibility (VIS 
Observations by VIS Code 


Category 


Category 


Weather 

(ww) 


Weather 

(W) 


Type 

(CL) 


All 

(90-99) 


Fair 

(94-95) 


Poor 

(90-93) 


Strong 


al 


41G 


4,5 


6 


66 


19 


32 


Foggers 


a2 


1 1 


If 


B 


466 


83 


360 


= S 


bl 


tl 


2G 


6 


24 


7 


2 




b2 


1! 


1 1 


B 


28 


6 


14 




b3 


10G 


4,5 


6 


51 


9 


2 




b4 


f I 


1 1 


B 


36 


12 


16 




cl 


41G 


* 


6 


2 


1 


1 




c2 


tl 


4,5 


5,7 


26 


5 


13 




c3 


tt 


It 


* 


64 


14 


39 




c4 


10G 


2G 


6 


30 


8 


2 


Total S cases: 








793 






Foggers 


di 


41G 


•JU 

A 


B 


21 


3 


17 


= F 


d2 


rt 


2G 


5,7 


12 


2 


3 




d3 


50G 


4 


6 


11 


3 


0 




d4 


1 1 


ft 


B 


13 


5 


6 




d5 


10G 


«ju 

A 


6 


3 


2 


0 




d6 


tt 


2G 


B 


13 


4 


1 




d7 


tt 


4 


5,7 


14 


2 


1 




el 


41G 


JL. 

A 


I f 


5 


2 


2 




e2 


ft 


2G 


* 


14 


4 


5 




e3 


50G 


4 


5,7 


3 


0 


0 




e4 


10G 


* 


B 


1 


0 


1 




e5 


tt 


4 


* 


37 


9 


1 




e6 


11G 


4 


6 


6 


0 


0 




el 


tt 


1 1 


B 


7 


2 


1 
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TABLE IV. (continued) 



fl 


41G 


* 


* 


2 


0 


1 


f2 


50G 


4 


J- 

/\ 


1 


0 


0 


f3 


10G 


5G 


5,7 


13 


3 


0 


f4 


11G 


ft 


6 


7 


0 


1 


f5 


It 


t f 


B 


12 


1 


0 


F gl 


10G 


Vc 


5,7 


2 


0 


1 


g2 


» » 


5G 


* 


24 


3 


0 


g3 


11G 


* 


6 


1 


0 


0 


g^ 


ft 


4 


5,7 


6 


2 


1 


g5 


ft 


f f 


J- 

A 


1 


1 


0 


hi 


ft 


J- 


B 


0 


0 


0 


h2 


f f 


5G 


5,7 


16 


2 


1 


il 


10G 


* 


* 


2 


0 


1 


i2 


11G 


JU 


5,7 


2 


2 


0 


jl 


ft 


5G 


* 


3 


2 


0 


Total F cases: 








262 






Past /Weak kl 


50G 


5G 


6 


61 


- 


8 


Foggers k2 


60G 


4 


6 


3 


- 


0 


II 

•"d 


f f 


ft 


B 


4 


- 


2 


£2 


* 


ft 


6 


32 


- 


2 


rnl 


50G 


' C* 


6 


6 


- 


1 


m2 


* 


4 


B 


31 


- 


12 


nl 


* 


5G 


6 


313 


- 


2 


ol 


60G 


4 


5,7 


1 


- 


0 


o2 


* 


4 


5,7 


39 


- 


0 


Pi 


60G 


4 


/> 


1 


- 


0 


p2 


JU 

/\ 


JU 


6 


26 


- 


0 


P3 


JU 


4 


JU 


61 


- 


1 


Total P cases: 








578 
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TABLE IV. (continued) 



Maybe ql 


50G 


5G 


B 


51 


- 


- 


Foggers q2 


f f 


tf 


5,7 


29 


- 


- 


= M q3 


ft 


* 


5,7 


7 


- 


- 


q4 


f f 


5G 


* 


19 


- 


- 


q5 


11G 


* 


k 


0 


- 


- 


q6 


60G 


5G 


6 


14 


- 


- 


q7 


11 


* 


6 


21 


- 


- 


rl 


50G 




B 


7 


- 


- 


r2 


tf 


/\ 


* 


3 


- 


- 


r3 


60G 


5G 


B 


5 


- 


- 


r4 


f ! 


tf 


5,7 


19 


- 


- 


r5 


ft 


A 


B 


19 


- 


- 


r6 


f ! 


* 


5,7 


50 


- 


- 


Total M cases: 








244 






Non-Foggers ul 


60G 


5G 


* 


8 


- 


0 


= N u2 


t f 


* 


* 


17 


- 


1 


vl 


* 


5G 


B 


192 


- 


2 


wl 


* 




B 


13 


- 


6 


xl 


A 




5,7 


44 


- 


9 


yi 


JU 

A 


5G 


5,7 


687 


3 


2 


y2 


JL 

A 


5G 


5,7 


1541 


5 


2 


zl 


* 


* 


* 


102 


1 


1 


Total N cases: 








2604 







Note: A dash in a visibility scheme indicates that particular count 

\jas not done. 
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TABLE V. Major fog categories and sub-category fog groups 
with associated Base and Fog Scores as a function 
of visibility. 

Visibility Code Values 





96-99 

(Good) 


94-95 

(Fair) 


90-93 

(Poor) 


Base 

Score 


Fog 

Score 


Maj or 
Category 


Subcategory Groups 












a 


580 


100.0 






a 


b 


560 


96.6 




a 


b 


c 


540 


93.1 


S 


b 


c 


d 


520 


89.7 




c 


d 


e 


500 


86.2 




d 


e 


f 


480 


82.8 




e 


f 


g 


460 


79.3 




f 


g 


h 


440 


75.9 


F 


g 


h 


i 


420 


72.4 




h 


i 




400 


69.0 




i 


i 




380 


65.5 




,1 




k 


360 


62.1 




k 


k 


J l 


340 


58.6 




l 


l 


m 


320 


55.2 




ra 


m 


n 


300 


51.7 


P 


n 


n 


o 


280 


48.3 




o 


o 


P 


260 


44.8 




P 






240 


41.4 





q 


q 


q 


220 


37.9 


M 


r 


r 


r 


200 


34.5 


not 
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TABLE VI. Correlation coefficients between predictor parameters and Fog Score (predictand) , 
and range and mean values of parameters, July 1976, North Pacific Ocean. 
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**Data base for these are the 4481 values interpolated to the reporting ship positions. 



TABLE VII. Regression equation predictor parameters and coefficients, and associated explained 
predictand variances (R 2 ) for stepwise selection, whole area, North Pacific Ocean 
(about 30-60 N) , July 1976. 
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TABLE VIII. Regression equation predictor parameters and coefficients, 
and associated explained predictand variances (R 2 ) for 
special treatment of climatology parameter (final-step 
form) , whole area, North Pacific Ocean (about 30-60 N) , 
July 1976. 

Equation 9T Equation 9NC 



Fog Score transformed (T) 

(Fog Score = Fog Score - CLUD) 
before regression. 



Regression 

Coefficient 


Variable 
(in entry- 
order) 


AR 2 


1 . 

-1.8328 


CUM) 

EHF 


0.0299 


+2.9127 


EAIR* 


0.0279 


-2.0207 


TAIR* 


0.0154 


40.5283 


V 


0.0046 


-1.3244 


ASTDR 


0.0022 


- 


YI upon entry 


0.0557 


YI upon removal 


-0.0004 


Final R 2 = 


0.1353 





FTER, DEWW, rSEHF, THF, SHF, 
CAPU, CAFV)** not in the use 
list. 



Regression 

Coefficient 


Variable 
(in entry 

order) 

* 

YI 


AR 


55.9906 


0.4420 


-1.8859 


•k 

EHF 


0.1183 


-4.4919 


ASTDR 


0.0035 


+21.0621 


ADTAIR 


0.0042 


-4.9522 


TSEA 


0.0028 


+2. 9131 


■k 

EAIR 


0.0091 


40.9988 


T9 25 


0.0028 


-0.2333 


PIW 


0.0012 


40.2698 


V 


0.0004 



Final R 2 = 0.5843 

CUM), FTER, DEWW, (SEHF, THF, 
SHF, CAPU, CAPV)** not in the 
use list. 



2 

* Most Significant Variables (based on AR > 0.005). 

** Omitted through error. 
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Regression equation predictor parameters and coefficients, 
and associated explained predictand variances (R 2 ) , final- 
step form, for two sub-areas, North Pacific Ocean, July 1976. 

Equation 40 

Based on 2977 observations at and 
north of 40 N. 



TABLE IX. 



Equation 39 

Based on 1504 observations south 
of 40 N. 



Regression 


Variable 


AR 2 


Coefficient (in entry 
order) 




-0.9826 


EHF^ 


0.0929 


0.6595 


CLIMO* 


0.0246 


0.5560 


EAIR* 


0.0089 


-1.4355 


SSTA 


0.0025 


-0.1608 


SEW 


0.0016 


-74.0931 


stabr* 


0.0018 


YI 


upon entry 


0.2842 


YI 


upon removal 


-0.0003 


HW 


upon entry 


0.0043 


HW 


upon removal 


-0.0014 


Final R 2 = 


0.4219 





Regression Variable 

Coefficient (in entry 
order) 


AR 2 


21.8849 


* 

YI 


0.5337 


-2.0569 


•k 

EHF 


0.0601 


-7.3865 


ASTDR* 


0.0058 


0.2945 


CLIMO 


0.0020 


-7.6194 


TSEA* 


0.0165 


0.7898 


T925 


0.0014 


0.5861 


V 


0.0011 


-26.5505 


AASTDX 


0.0010 




ADTSEA 
upon entry 


0.0047 




ADTSEA 
upon removal 


-0.0004 


Final R 2 = 


0.6278 





FIER,DDWW, J^SEHF, THF, SHF, 

CAPU, CAFV) not in use list 
for both equations. 

2 

* Mast Significant Variables (based on AR > 0.005). 

** Omitted through error. 

@ Not best equation (39) possible. A full use list run (not shown) 
made late in the study gave SEHF a higher initial correlation to 
Fog Score. 
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TABLE X. 


Heidke Skill and Threat Scores, whole area, North 
Pacific Ocean (about 30-60 N) , July 1976. FOGCAT 
designators S and F used to identify marine fog 
observations. See Table IV and text. 


Predictor 


Cutoff 

(Cut) 

(Per cent) 


Skill Score 


Threat Score 


FTER 


30 


0.274 






0.304 






70 


.029 






.023 


CLIMO 


30 


.271 






.325 






40 


.232 






.241 


Eq 


1 


40 


.340 






.355 


Eq 


2 


40 


.355 






.372 


Eq 


3 


40 


.346 






.365 


Eq 


4 


40 


.350 






.367 


Eq 


5 


40 


.361 






.372 


Eq 


6 


40 


.376 






.388 


Eq 


7 


40 


.378 






.389 


Eq 


8* 


40 


.386 






.394 


Eq 


9 


40 


.381 






.392 






45 


.383 






.373 


Eq 


9.9 


35 


.380 






.372 


Eq 


9T 


45 


.383 






.373 


Eq 


9NC 


40 


.369 






.384 



Scoring performed on whole area data file; i.e., 4481 cases, 
of which 1055 are classified Foggers, using major categories 
S and F as criteria. 

Error in one variable's coefficient makes these results 
slightly suspect. 



TABLE XI. Daily Heidke Skill Scores, whole area, North 
Pacific Ocean (about 30 - 60 N) , July 1976. 
FOGCAT designators S and F used to identify 
marine fog observations. See Table IV and text. 



DAY 


Cases 


Foggers 


FTER 
Cut=30 
(Per cent) 


CLIMO 
= 30 


Eq. ( ( 
= 40 


1 


125 


28 


0.159 


0.336 


0.459 


2 


162 


34 


.119 


.346 


.273 


3 


145 


35 


.140 


.245 


.282 


4 


147 


33 


.132 


.228 


.379 


5 


152 


39 


.016 


.149 


.214 


6 


156 


30 


-0.060 


.200 


.172 


7 


148 


27 


.027 


.160 


.331 


8 


140 


33 


.050 


.322 


.332 


9 


160 


32 


.219 


.326 


.408 


10 


163 


38 


.459 


.372 


.540 


11 


132 


25 


.284 


.285 


.315 


12 


145 


33 


.319 


.331 


.442 


13 


159 


32 


.231 


.190 


.213 


14 


136 


21 


.010 


.172 


.211 


15 


160 


34 


.323 


.386 


.455 


16 


156 


34 


.228 


.328 


.363 


17 


141 


31 


.238 


.200 


.336 


18 


164 


38 


.292 


.148 


.402 


19 


165 


40 


.414 


.264 


.513 


20 


147 


32 


.450 


.303 


.446 


21 


138 


24 


.118 


.162 


.291 


22 


151 


40 


.407 


.416 


.438 


23 


148 


40 


.428 


.378 


.483 


24 


137 


32 


.488 


.386 


.539 


25 


143 


45 


.444 


.367 


.512 


26 


156 


42 


.455 


.490 


.451 


27 


146 


41 


.342 


.251 


.337 


28 


143 


44 


.348 


.107 


.337 


29 


157 


49 


.505 


.212 


.402 


30 

Total 


159 

4481 


49 

1055 


.267 


.161 


.421 



CLIMO has more skill on 17 of the 30 days compared to FTER. 
Eq.(9) has more skill on 24 of the 30 days compared to FTER. 

Eq.(9) has more skill on 28 of the 30 days compared to CLIMO 

Eq.(9) has more skill on 23 of the 30 days compared to CLIMO 

and FTER 
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TABLE XII. Daily Threat Scores, whole area, North Pacific 
Ocean (about 30-60 N) , July 1976. FOGCAT 
designators S and F used to identify marine fog 
observations. See Table IV and text. 



DAY 


FTER 




CLIMO 


EQ. (9) 








Cut=30 




=30 


=40 








(Per cent) 










1 


0.241 




0.346 


0.446 






2 


.202 




.301 


.306 






3 


.247 




.321 


.342 






4 


.229 




.299 


.385 






5 


.198 




.276 


.314 






6 


.134 




.257 


.247 






7 


.136 




.224 


.328 






8 


.194 




.349 


.350 






9 


.231 




.328 


.374 






10 


.421 




.377 


.491 






11 


.283 




.305 


.321 






12 


.339 




.349 


.421 






13 


.229 




.250 


.241 






14 


.135 




.219 


.231 






15 


.316 




.379 


.431 






16 


.262 




.355 


.375 






17 


.233 




.269 


.362 






18 


.265 




.250 


.400 






19 


.354 




.325 


.485 






20 


.412 




.346 


.424 






21 


.170 




.229 


.300 






22 


.400 




.426 


.439 






23 


.429 




.408 


.477 






24 


.464 




.293 


.509 






25 


.471 




.427 


.539 






26 


.457 




.486 


.458 






27 


.407 




.353 


.414 






28 


.388 




.284 


.410 






29 


.514 




.347 


.451 






30 


.314 




.304 


.456 






CLIMO has more 


skill on 19 


of 


the 30 days 


compared 


to 


FTER. 


Eq. (9)has more 


skill on 29 


of 


the 30 days 


compared 


to 


FTER. 


Eq . (9) has more 


skill on 27 


of 


the 30 days 


compared 


to 


CLIMO 


Eq . (9) has more 


skill on 26 


of 


the 30 days 


compared 


to 


CLIMO 










and 


FTER. 
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TABLE XIII 



. Heidke Skill and Threat Scores, sub-areas, North 
Pacific Ocean (about 30-60 N) , July 1976. FOGCAT 
designators S and F used to identify marine fog 
observations. See Table IV and text. 



Predictor 
Area : 


Cut 

(Per cen£ 
Northern Belt 


Skill Score 

) 

(2977 observations; 


Threat Score 
948 foggers) 


FTER 


30 


0.233 


0.335 




70 


.026 


.026 


CLIMO 


30 


.151 


.341 




40 


.170 


.259 


Eq. (9) 


40 


.299 


.412 




45 


.322 


.397 


Eq. (40) 


40 


.306 


.418 




45 


.337 


.412 



Area: 


Southern Belt' 


~(1504 observations; 


107 foggers) 


FTER 


30 


0.156 


0.133 




70 


-0.001 


.000 


CLIMO 


30 


.067 


.043 




40 


.000 


.000 


Eq. (9) 


35 


.288 


.213 




40 


.247 


.174 


Eq. (39) 


30 


.283 


.209 




35 


.254 


.173 



* Northern Belt = whole area cases at or north of 40 N 
^Southern Belt = whole area cases south of 40 N 
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TABLE XIV. 



Probability Scores (PS) for whole and sub-areas 
of the North Pacific Ocean (30-60 N) , July 1976. 
FOGCAT designators S and F used to identify 
marine fog observations. See Table IV and text. 



Whole Area 



Predictor 



PS 



Max 



Min 



No. of observations 
with value below zero 



FTER 



CLIMO 

Eq. (1) 
Eq. (2) 
Eq. (3) 
Eq. (4) 
Eq. (5) 
Eq. (6) 
Eq. (7) 
Eq. (8)* 
Eq. (9) 
Eq. (9.9) 
Eq. (9T) 
Eq. (9NC) 



0.350 


1.02 


-0.09 


1183 


)2 .349 


1.04 


-0.07 


292 


)5 .347 


1.07 


-0.04 


21 


.315 


59.3 


0.0 


0 


.323 


96. 


-48. 


195 


.319 


87. 


-36. 


195 


.317 


87. 


-33. 


185 


.317 


88. 


-35. 


209 


.315 


94. 


-37. 


192 


.313 


85. 


-53. 


257 


.312 


88 


-53. 


262 


.312 


88. 


-52. 


268 


.311 


89. 


-50. 


253 


.290 


79. 


-60. 


576 


.313 


88. 


-55. 


261 


.314 


87. 


-47. 


277 


Belt (cases 


at or north of 40 N) 




Predictor 


PS 


Predictor 


PS 


FTER 


0.446 


Eq. (9) 


0.393 


CLIMO 


0.412 


Eq. (40) 


0.388 


Belt (cases 


south of 


40 N) 




FTER 


0.162 


Eq. (9) 


0.149 


CLIMO 


0.124 


Eq. (39) 


0.135 



c 

Error in coefficient makes result slightly suspect, 
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TABLE XV. Daily Probability Scores, whole area, North Pacific 
Ocean (about 30-60 N) , July 1976. FOGCAT designa- 
tors S and F used to identify marine fog observa- 



tions. See 



DAY 




FTER 


1 




0.384 


2 




.391 


3 




.394 


4 




.434 


5 




.471 


6 




.447 


7 




.350 


8 




.388 


9 




.341 


10 




.281 


11 




.299 


12 




.340 


13 




.302 


14 




.305 


15 




.310 


16 




.355 


17 




.366 


18 




.365 


19 




.347 


20 




.364 


21 




.312 


22 




.334 


23 




.337 


24 




.267 


25 




.311 


26 




.307 


27 




.332 


28 




.369 


29 




.344 


30 




.450 


CLIMO 


had 


more skill 


Eq. (9) 


had 


more skill 


Eq. (9) 


had 


more skill 


Eq. (9) 


had 


more skill 


Eq. (9.9)had 


more skill 


Eq. (9.9)had 


more skill 


Eq. (9 . 9)had 


more skill 


Eq. (9.9)had 


more skill 



Table IV and text. 



CLIMO 


Eq. (« 


0.281 


0.291 


.318 


.317 


.341 


.352 


.305 


.328 


.369 


.385 


.299 


.342 


.279 


.318 


.314 


.332 


.276 


.285 


.301 


.284 


.270 


.303 


.291 


.298 


.309 


.300 


.259 


.282 


.265 


.271 


.299 


.316 


.308 


.320 


.328 


.298 


.314 


.276 


.305 


.278 


.282 


.299 


.337 


.295 


.326 


.292 


.370 


.252 


.351 


.308 


.297 


.291 


.366 


.368 


.390 


.362 


.394 


.335 


.391 


.348 



on 22 of the 30 days 
on 27 of the 30 days 
on 14 of the 30 days 
on 13 of the 30 days 

on 29 of the 30 days 
on 30 of the 30 days 
on 27 of the 30 days 
on 26 of the 30 days 



Eq. (9.9) 

0.272 

.289 

.326 

.293 

.364 

.246 

.272 

.308 

.260 

.271 

.262 

.270 

.280 

.245 

.253 

.288 

.296 

.288 

.265 

.259 

.257 

.297 

.289 

.228 

.307 

.282 

.341 

.355 

.336 

.354 



compared to FTER. 
compared to FTER. 
compared to CLIMO. 
compared to CLIMO 
and FTER. 
compared to FTER. 
compared to CLIMO. 
compared to Eq . (9) 
compared to FTER, 
CLIMO, and Eq . (9) 



TABLE XVI. Total observation counts, and fog observation counts and frequencies, 
by probability class, whole area (about 30-60 N) , July 1976. FOGCAT 
designators S and F used to identify marine fog observations. See 
Table IV and text. 
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TOTAL 4481 1055 23.5 4481 1055 23.5 4481 1055 23 







Figure 1. Study Area on Mercator projection. Land stations used in 
study and Ocean-Weather-Ship P are encircled. 
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Figure 2. Study Area on Polar Stereographic projection. Grid point 
locations involved in advection computations requiring 
bogusing are labeled. See Appendix D. 




Figure 3. Fleet Numerical Weather Central's 63x63 grid 
with outline of North Pacific Ocean rectangu 
lar grid area used in study. See text. 
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Figure 4. Naval Postgraduate School's fog frequency climatology (percent) 
for July. (After Willms, 1975). 
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Figure 5. National Meteorological Center's sea-level pressure analysis 
(tens, units, mb), 0000 GMT, 12 July 1976. 
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Figure 6. Analysis of FNWC's evaporative heat flux grid point data 
(g • °C/cm'- • depth) for 0000 GMT, 12 July 1976. 
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Figure 7. Analysis of FNWC ' s fog probability (FTER) grid point 'data 
(per cent) for 0000 GMT, 12 July 1976. 
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Figure 8. Analysis of Fog Score grid point data [from regression 
equation (9)] (per cent) for 0000 GMT, 12 July 1976. 
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Figure 9. Categorized synoptic observations (FOGCAT) for 0000 GMT, 
12 July 1976. S = Strong Foggers, F = Foggers , P = Past/ 
Weak Foggers, M = Maybe Foggers, N = Non-Foggers. 




83 



Figure 10. National Meteorological Center's sea-level pressure analysis 
(tens, units, mb), 0000 GMT, 24 July 1976. 
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Figure 11. Analysis of FNWC's evaporative heat flux grid point data 
(g • °C/cm2 • depth) for 0000 GMT, 24 July 1976. 
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Figure 14. Categorized synoptic observations (FOGCAT) for 0000 GMT, 
24 July 1976. S = Strong Foggers , F = Foggers, P = Past/ 
Weak Foggers, M = Maybe Foggers, N = Non-Foggers. 
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